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We use recently developed effective field theory techniques to calculate the third order post- 
Newtonian correction to the spin-spin potential between two spinning objects. This correction 
represents the first contribution to the spin-spin interaction due to the non-linear nature of general 
relativity and will play an important role in forthcoming gravity wave experiments. 
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The problem of motion in general relativity has shed 
its seeming academic nature due to pending gravity wave 
experiments. Both LIGO [jj] and LISA [^ expect to de- 
tect radiation from inspiralling binaries, and thus build- 
ing templates for these events has become increasingly 
important. While for late stages of the inspirals numer- 
ical techniques are needed, for the early stages one may 
calculate in the Post-Newtonian approximation (PN) for 
small velocities. In relativistic cases one may also calcu- 
late analytically when deviations from the Schwarzschild 
geometry are small. A complicating factor in the build- 
ing of these templates is the fact that there are mul- 
tiple scales involved. In particular, the finite size of 
the object, leads to tidal deformations and dissipation 
which can then in turn affect potentials as well as ra- 
diation. These effects make an exact analytical solution 
intractable. One might hope that systematically expand- 
ing around the point particle approximation would lead 
to a controlled calculational scheme. This is in fact what 
one can do by using techniques developed for effective 
field theories |j, |j, |^. While these theories are typi- 
cally applied to quantum field theories involving multiple 
scales, we will be applying them here to a purely classical 
problem. The classical problem shares many of the same 
calculational hurdles as the quantum problem. In partic- 
ular, because we are expanding around the point particle 
approximation, the calculation of potentials, as well as 
the radiated power loss, entails regularizing divergences. 
However, since in the EFT we work at the level of the 
action, these divergences are simply renormalized|3, and 
for certain higher dimensional terms in the action lead 
to a classical renormalization group trajectory jg. The 
EFT also allows one to power count in a very natural 
way. Each term in the action has a definite scaling in 
the expansion parameter and thus we may determine the 
size of a given effect simply by reading off the scaling of 
the relevant terms in the action. 

In the last few years the 3PN {{0{v^)) potentials for 
non-spinning objects have been computed 13, but the 
case of spinning bodies has not reached that level of ac- 
curacy. The 2.5PN spin-orbit potential was calculated in 
[8| within a point particle approximation and Hadamard 
regularization. The leading order spin potentials were 
obtained in |2|, but the spin-spin piece has yet to be 
computed to 3PN. By using the techniques developed in 



[2j we may avoid any of the pitfalls of the point parti- 
cle approximation and tame divergences as well as finite 
size effects. Following the literature, in the spinning case 
when we use the term 3PN, we mean order v^v'^^^, where 
Vrot is the rotational velocity which is taken to be or- 
der one, in the "maximally rotating case" and ev in the 



"co-rotating case", where e = R/r ('■ 



for compact 



objects), and R, r are the radius of the object and orbit 
respectively. 

In this letter we will be appl ying these EFT ideas, ex- 
tended to spinning particles in 10], to calculate the first- 
nonlinear correction to the spin-spin potential for com- 
pact binaries which arises at 3PN ^. This correction is 
the "hyper-fine" analog of the Einstein-Infeld-Hoffmann 
potential calculated more than seventy years ago 11| . 
While the high order of this contribution makes it ap- 
pear to be of only academic interest, such accuracy may 
indeed be necessary for successful matching 12] due to 
the fact that inspirals are tracked over long periods of 
time. 

As explained in [2j , the EFT approach proceeds in two 
stages. First one matches the full theory of an extended 
object interacting with gravity onto a world-line action 
treating each object in isolation. This action consists 
of a tower of world-line operators, with the coefficients 
of higher dimensional operators encapsulating the true 
finite size nature of the underlying particle. These coeffi- 
cients can be fixed by a matching procedure as discussed 
for the case of absorption in (J|. At the order we are 
working the world line action necessary for our calcula- 
tion is given by 

S= - Y. f^adSa+ f ^Syn'}j{Xa)dXa 

- Ml I d^x [2^gR + {d^h^, ~ \d^hlf^ (1) 

where Sa and Xa are the proper length and affine param- 
eters for the a'th worldline respectively. The last term is 
the gauge fixing term, which corresponds to the harmonic 



^ The results can be easily extended to non-compact objects. In 
the latter the leading and subleading contributions to the spin- 
spin potential scale as tv^^^ and « tij;^^!) respectively. 



gauge used in all our calculations. In llOl . following the 
classic work of Regge and Hanson [13 for a relativistic 
top in flat space, the spin effects are included in a gener- 
ally covariant fashion by introducing the vierbein degrees 
of freedom e^{\a) on the worldline, which relate the local 
co-rotating frame to the global frame. The generalized 
angular velocity is then given by (for each particle) 






(2) 



and the spin Sjj is the variable conjugate to Qij. Op- 
erators describing finite size effects have been left out as 
they are subleading corrections to the spin-spin effects 
we are interested here [j, |li| . We will come back to this 
point later on. The form of the world-line action (^ is 
fixed by reparametrization invariance, which allows us to 
freely choose the worldline parameter. A wise choice is 
thus X = x^, which will directly lead us to the effective 
action in the PN frame. Therefore, throughout this letter 

We next match onto a theory of potentials. This is 
accomplished by first expanding the action around the 
flat space limit, i.e 






= Vt^y + h 
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2M, 



pi 






(3) 



where A is the Lorentz transformation that relates the 
local co-rotating and global frames in the flat space limit. 
Using the power counting rules developed in [3, llOf we 
expand each term in powers of the relative velocity such 
that each term in the action scales homogeneously in v. 
This allows us to calculate potentials to arbitrary order 
in V by drawing all possible Feynman diagrams involving 
a flxed set of vertices. Before specifying the action, we 
must choose a gauge for the spin degrees of freedom as 
well as the tetrad. A common and convenient choice are 
the Newton- Wigner coordinates 



60 



i^°'^ 



m 



S-'O = -v^S'^ +0{v^) 



(4) 



in which the position and momenta are treated as canon- 
ical variables and the spin obeys the traditional angular 
momentum algebra (l^. These constraints, which are not 
all linearly independent, eliminate the unphysical degrees 
of freedom. Results in this gauge can be related to other 
choices via the appropriate coordinate transformations 

Given this gauge choice, we expand out the action 
yielding to the leading order spin graviton coupling 



L2 = 



1 



2Af, 



-S'^h°''\ 



pi 



(5) 



which arises from expanding fi in terms of the tetrad and 
using ^. Given that the spin scale as [lOJ 



5' ex 



m 



M, 



-Vrot, 



pP 



(6) 



we find that this leading order interaction goes like v^Vrot- 
Thus the leading order spin-spin potential follows from 
the Feynman diagram with two of these leading order 



interactions and thus will scale as 



It is given by 



a one graviton exchange diagram leading to 



^2PN — 



— 5-(<3i 



^2 



J ■ Sir ■ S2 



(7) 



which agrees with the previously derived results 9] . We 
will use throughout the paper 5'^ = e^^^S^ (for differ- 
ent definitions see [8|, ll3j)j and the standard notation 
labelling the post-Newtonian order by the power of the 
orbital velocity. 

To work to next order we must consider vertices in- 
volving spin to order v^Vrot- At this order we have 



L.= 



1 



2Af, 



.^.jfe^^/j^fcj 



1 



pi 



AM, 



-v'^^'/i™^^ 



pi 



(8) 



while at order v'^Vrot we must include 



L. 



1 



-v^S^^v'h'°'^ + 



1 



-v'S'^h°^'°. 



AMpi AMpi 

These terms arise by keeping higher order pieces in the 
spin matrix (i.e. Soi) as well as time derivatives which 
are down by v relative to spatial derivatives. Since we 
are calculating to order v^v^^f (3PN) we must include 
diagrams with one insertion of L4 and a leading order 
vertex L2, as well as those diagrams with two insertions 
of L3. These diagrams are depicted in Figures (la) and 
(lb). There is a further contribution from one graviton 
exchange which arises from the first correction to the 
graviton Green's function 



G(p) ex -2- 
Po 



1 



jp 



1 

p2 



Po 

P" 



(9) 



which reflects the first deviation from instantaneity in 
the exchange. This contribution is depicted in Figure 
(Ic) and^ more formally, arise from new vertices in the 
action '3\. 

At order 3PN we must also include diagrams with dou- 
ble graviton exchange (note that we need not include di- 
agrams which can be disconnected by cutting the world- 
line, as these are iterations of connected diagrams which 
are automatically resummed when solving the equations 
of motion) . Terms quadratic in the metric arise from the 
higher order terms in the vierbein and the connection. 
After expanding we find that the relevant pieces of the 
Lagrangian are 



1 



^/l2 



^Kl 



S'^{h'''{h^^'^ - W^^^) - h'%"°'^). (10) 




b) 



c) 



FIG. 1: Diagrams contributing to 3PN order which do not 
involve non- linear interactions. The blow represents a spin in- 
sertion and the cross corresponds to a propagator correction. 
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FIG. 2: Non-linear contributions to the 3PN spin-spin poten- 
tial. 



These terms scale as v'^Vrot thus, we must include all 
diagrams for which a leading order spin vertex {v^Vrot) 
and a leading order mass insertion from the non-spinning 
part of the action 



io = - 



2M 



pi 



(11) 



are contracted. The resulting diagram is shown in Fig- 
ure (2b). Finally we must include the contribution from 
those diagrams including the three graviton vertex which 
scales as w^ |3] . This vertex is too voluminous to include 
here, and is best, and quite simply, handled using a sym- 
bolic manipulation program All such diagrams resulting 
from this vertex, two leading order spin insertions and 
one leading order mass insertion must be included (see 
Figure (2a)). 

Combining all of the contributions leads to the poten- 
tial 



V, 



_ Gn 

3PN - -^ 



Si ■ S2 { -Vi ■ V2 - 3ui • TIV2 ■ n 



Si ■ V1S2 ■ V2 — -;;Si ■ V2S2 ■ vi + Si ■ V2S2 ■ V2 



+ S2 ■ viSi ■ vi + 35i • nS2 ■ n{vi ■ V2 + 5vi ■ nv2 ■ ft) — 35i ■ V1S2 ■ nv2 ■ n — 35*2 • V2S1 ■ nvi ■ n 

^ 3 ^ 

-I- 3(u2 X Si) ■ n{v2 X S2) ■ n + 3(ui x Si) ■ n{vi x 6*2) • n - -{vi x Si) ■ n{v2 x S2) ■ n 



6(wi X ^2) • n{v2 X 5*1) • n 



3G%{Tni + 1712) 



Si ■ S2 — 3Si ■ nS2 ■ n 



(12) 



with ft = -. The equations of motion simply follow from 
the usual Hamiltonian procedure. 

The leading order finite size contribution arise from so- 
called self induced effects. That is, corrections to spheric- 
ity arising from the non-vanishing quadrupole moments 
induced by rotation. These effects are encapsulated by 
the world-line operator 



C 



^ES^ 



2mM, 



'ipi 



2 P 



(13) 



where C is a constant which is determined by the nature 
of the object and E^i, is the electric component of the 
Wcyl tensor. In the case of a rotating black hole C = 1, 
and this term represents the non-vanishing quadrupole 
moment of the Kerr solution [l5j. For neutron stars C 
ranges between 4 and 8 depending on the equation of 
state of the neutron star matter [la| • This operator scales 
as v'^v^gf at leading order, and it is easy to show that 
it gives rise to a 2PN (quadrupole-monopole) spin-orbit 



contribution pl(ill7|. 



CGNni2 
2mir^ 



^2PN — 



Si ■ Si - 3{n ■ Si)^) +1^2, 



(14) 

so one might then think that it will also contribute at 
3PN in the spin-spin sector, since naively one could con- 
tract H13|l with a leading order spin interaction |(SJ| . How- 
ever, this contribution vanishes, and the leading order 
finite size effects in the spin-spin potential show up at 
3.5PN. There is nevertheless a finite size correction in 
the spin-orbit interaction at 3PN coming from diagrams 
where (jl3|l is contracted with subleading mass insertions 
as well as corrections to instantaneity and non linear ef- 
fects. We will report on these results in a forthcoming 
paper. 

As discussed in J3| the inclusion of radiation into this 
formalism is accomplished by working in a background 
field. In the background field gauge, after multipole ex- 
panding the field to keep manifest power counting |18l |. 
and integrating over (i.e. solving the equations of mo- 
tion) the potential gravitons, one generates an effective 



action for the sources and the radiatfon graviton r[a;, i, /i] 
which is invariant under small diffeomorphisms. The 
power loss due to radiation then follows by calculating 
the imaginary parts of all vacuum diagrams (i.e. no ex- 
ternal gravitons), whereas the real part introduces ra- 
diation reaction effects into the equations of motion for 
the sources. Including spin effects in radiation follows 
in a similarly straightforward manner by simply includ- 
ing vertices in which the spin couples to the background 
field. We have calculated the effects of spin on radiation 
to 2PN and find agreement with those in the literature 
[19j. These results along with new results for the 3PN 
contributions to radiation and quasicircular orbits with 
the inclusion of spin, will follow in a subsequent publica- 
tion. 

Finally, it is hopefully clear that there are no stum- 
bling blocks in going beyond 3PN in our formalism. It 
is simply a matter of bookkeeping. Since every term 
in the action scales homogeneously in the velocity this 
is not a difficult task. One simply draws all possible 
diagrams such that the net scaling of all the vertices is 
of the desired accuracy. For instance, at next order one 
would need the four-graviton vertex as well as the higher 
order vertices with and without spin and a contact 
term involving one source and three graviton lines. 
Divergences are not an issue as they are all absorbed 
into their appropriate counter-terms. They are irrelevant 
until one reaches the order at which finite size operators 
become relevant, at which point rcnormalization is 
straightforward. 
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